Kovacs effects in an aging molecular liquid 
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We study by means of molecular dynamics simulations the aging behavior of a molecular model 
of ortho-terphenyl. We find evidence of a a non-monotonic evolution of the volume during an 
isothermal-isobaric equilibration process, a phenomenon known in polymeric systems as cross-over 
(or Kovacs) effect. We characterize this phenomenology in terms of landscape properties, providing 
evidence that, far from equilibrium, the system explores region of the potential energy landscape 
distinct from the one explored in thermal equilibrium. We discuss the relevance of our findings for 
the present understanding of the thermodynamics of the glass state. 
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If two systems in thermodynamical equilibrium with 
identical chemical composition have the same tempera- 
ture T and volume V we immediately know that they 
experience the same pressure P. We know that the two 
systems will respond in the same way to an external per- 
turbation, and that they will be characterized by the 
same structural and dynamical properties. The ability 
to predict the pressure, and the equivalence of structural 
and dynamical properties, derive from thermodynamical 
principles. 

In the case of glasses, systems in out-of-equilibrium 
conditions, the T and V values are not sufficient for 
predicting P, since the state of the system depends on 
its previous thermal and mechanical history. Different 
glasses, at the same T and V, are characterized by differ- 
ent P values. One can ask if the value of P is sufficient 
to uniquely define the glass state, i.e., if two glasses with 
identical composition having not only the same T and V 
but also the same P are the same glass. If this is the 
case, the two glasses should respond to an external per- 
turbation in the same way and should age with a similar 
dynamics. 

These basic questions are at the hearth of a thermody- 
namic understanding of the glassy state of matter, and of 
the possibility of providing a theoretical understanding 
of out-of-equilibrium systems. Indeed, if by specifying 
T, V and also P, we uniquely define the glass state — its 
structural and dynamical properties — it means that it is 
possible to develop an out-of-equilibrium thermodynamic 
formalism 0, IE IE El IS El where the previous history of 
the system is encoded in one additional parameter. In 
the interpretation of experimental data, such additional 
parameter is often chosen as a fictive temperature or pres- 
sure, in the attempt to associate the glass to a liquid, 
frozen from a specific thermodynamic state. 

Back in the 60th, Kovacs and co-workers designed an 
experimental protocol 0, H, 13 (Fig. to generate dis- 
tinct glasses with different thermal and mechanical histo- 
ries but with the same T, V and P values. Poly-vinyl ac- 
etate was equilibrated at high temperature Th and then 
quenched at low temperature Tj, where it was allowed 
to relax isothermally for a waiting time t w insufficient 



to reach equilibrium. The material was then re-heated 
to an intermediate temperature T, and allowed to re- 
lax. The entire experiment was performed at constant 
pressure P. The observed dynamics of the volume relax- 
ation toward equilibrium — in the last step at constant 
T and P — was striking; the volume crosses over the 
equilibrium value, passes trough a maximum, which de- 
pends upon the actual thermal history of the system, and 
then relaxes to the equilibrium value. The existence of a 
maximum clearly indicates that there are states with the 
same V (at the left and at the right of the maximum) 
which, although T, V and P are the same, evolve differ- 
ently Thus, this experiment strongly support the idea 
that three variables are not always sufficient to uniquely 
predict the state of the glass [lOj . 

Here we attempt to reproduce numerically the Kovacs 
experiment, performing molecular dynamics simulations 
for a simple molecular model, to develop an intuition on 
the differences between states with the same T,V and P 
and the conditions under which out-of-equilibrium ther- 
modynamics may be used to describe glass states. We 
find that for sufficiently deep quenching temperatures, 
and long aging times, following the protocol proposed 
by Kovacs, it is indeed possible to identify two distinct 
states with the same T, V and P. Thus, for the first time, 
Kovacs' effects, also know in the literature as cross-over 
effects [n| , are also observed in a molecular liquid model. 
Exploiting the possibilities offered by the information en- 
coded in the numerical trajectories, we examine the dif- 
ferences between this pair of states, and contrast them 
with corresponding equilibrium liquid configurations. We 
discover that when the system is forced to age following 
significant T-jumps, i.e., low TJ, it starts to explore re- 
gions of the landscape which are never explored in equi- 
librium. Under these conditions, it is not possible any 
longer to associate a glass to a "frozen" liquid configu- 
ration via the introduction of a fictive temperature or 
pressure. This finding limits the range of validity of re- 
cent theories for out-of-equilibrium systems, based on the 
possibility of developing a thermodynamic formalism for 
glasses introducing only one additional effective parame- 
ter in the free energy 0, 0, IS IE 113 • 
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FIG. 1: Top: (a) Temperature protocol in Kovacs experiment : 
in the numerical simulation the system equilibrated at T;, is 
quenched at several Ti (1) where it relaxes for a time t w (2). It 
is finally heated at T (3). (b) Corresponding volume evolution 
in the cases of low (A) and high (B) Ti. Bottom: Volume 
relaxation at different Ti . The fast relaxation taking place for 
times shorter than the thermostat characteristic time is not 
reported. 



We consider a system of 343 molecules, interacting via 
the Lewis and Wahnstrom potential 13] (LW), a model 
for the fragile glass former ortho-terphenyl (OTP). The 
molecule is rigid, composed by three sites located at the 
vertexes of an isosceles triangle. Sites pertaining to dif- 
ferent molecules interact by the Lennard-Jones potential. 

Simulation details are given in Refs. mm. The sys- 
tem has been studied in isobaric-isothermal conditions; 
the time constants of both the thermostat and the baro- 
stat have been fixed to 20 ps. In order to accumulate an 
accurate statistics for the analysis below, averages over 
up to 200 different starting configurations have been per- 
formed. The total simulation time is more than 6 fis. 

To reproduce the Kovacs experiment, equilibrium con- 
figurations at Th = 400 K, volume per molecule V = 
0.378 nm 3 and P = 16 MPa are isobarically quenched 
at several low temperatures 7], and left to age for differ- 
ent t w values (t w is never sufficient to reach equilibrium 
at Ti). Each resulting configuration is then isobarically 
heated to T = 280 K (see Fig. [TJ. At T the charac- 
teristic structural relaxation time is of the order of a 
few ns, allowing us to follow, with the present compu- 
tational resources, the dynamics up to equilibrium. The 
1^-evolution is recorded along the entire path. We also 
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FIG. 2: Volume relaxation showing Kovacs effect in a molecu- 
lar system. Top: A system of 343 molecules of the LW OTP at 
constant pressure P — 16 MPa is equilibrated at Th = 400K, 
quenched at several low temperatures Ti and leaved to age 
for t w = 50 ns, a time insufficient to reach equilibrium at T; . 
The system is then heated at the intermediate temperature 
T = 280 K and the V relaxation dynamics is recorded (sym- 
bols). For the case T; = T, the system is directly brought 
from T h = 400A to T = 280A' and hence t w = 0. The verti- 
cal dashed lines indicate the transient time, where T and P 
have not yet equilibrated to the final values. Bottom: Volume 
relaxation at T — 280 K, at fixed Ti = 150 K for different t w 
values. Dashed lines delimitate the transient region, when 
neither T nor P have reached their equilibrium value. 



perform an analysis of the properties of the explored po- 
tential energy landscape (PEL) as a function of time, 
focusing in particular on the energy e/s and pressure 
Pis of the closest local minima configurations (the in- 
herent structures, IS), as well as the average curvatures 
of the PEL around the IS configurations. The IS con- 
figuration, which is numerically evaluated performing a 
steepest descent path along the potential energy surface, 
can be though of as the low-T glass generated by instan- 
taneously freezing the liquid. 

Fig. |2fa) shows the time evolution of V at constant 
T = 280 K and P = 16 MPa for samples which have 
previously aged at different T\ for t w = 50 ns. The time 
evolution of V for the case Ti — T (i.e. for a T-jump 
from Th to T) is also reported. We note that for large 
Ti (i.e., Ti = 220 K), V relaxes to the equilibrium value 



FIG. 3: Kovacs effect for T\ = 150 K and t w = 25 ns; aver- 
age over 200 independent realization. The out-of-equilibrium 
evolution of (a) system volume V, (b) inherent structures eis, 
(c) inherent structure pressure Pis, and (d) shape factor S 
are shown. Error bars are smaller than the symbols. Two 
points (labeled P < and P<), characterized by the same values 
of T, P, and V are marked by arrows. 



from below monotonically. Similarly, for TJ = T, V re- 
laxes to the equilibrium value from above, again mono- 
tonically. At variance, for deep Ti values, V does not 
relax monotonically to the equilibrium value. In anal- 
ogy with Kovacs findings, it goes through a maximum, 
whose value is higher the lower Ti, before starting to re- 
lax toward the equilibrium value. After the maximum, 
the time evolution of V practically coincides with case 
Ti=T. 

The waiting time dependence of V during the final re- 
laxation at T is shown in Fig. EJb) for the case T\ = 150 
K. For all studied t w values, a maximum is observed, 
and the value of the maximum is larger the shorter the 
waiting time t w . All together, Figs. |2{a) and (b) con- 
firm that the features observed by Kovacs in experiments 
with polymers are also observable in the case of molec- 
ular glass forming liquids under deep-quench conditions 
(small Ti values). 

Next we study the properties of the region of the PEL 
explored by the system during the the aging process. In 
Fig. |3 we show — for the case Ti = 150 K for which 
a clear non-monotonic V-relaxation is observed — the 
out-of-equilibrium evolution of (a) the V, (b) the aver- 
age inherent structures energy e/s, (c) the average in- 
herent structure pressure Pis, and (d) the shape factor 
S = J2k=i 3 fog^fc/^o)/^ (here the uju are the eigenval- 
ues of the Hessian calculated at the inherent structures 
and lo is the frequency unit); the quantity S provides a 
measure of the volume of the basin of attractions) [l6| . 
The maximum in the time evolution of V allows us to de- 
fine two arbitrary points P< and P> (marked in Fig. [2J , 
characterized, by construction, by the same T, P, and 
V values. Although P < and P> are indistinguishable 



-228 




0.352 



0.353 



0.354 



0.355 



V [ nm ] 



FIG. 4: Volume dependence of eis(a) ,Pis (b) and S (c) 
during the relaxation dynamics at T = 280K, for samples 
which have aged at T t = 150 K, 220 K and 280 K. Only data 
for times longer than 20 ps — the time associated to the es- 
tablishment of constant T and P conditions — are shown. 
Note that Pis is constant. The solid lines show the equi- 
librium relation eis(V), Pis(V) and S(V) along a constant 
Pis = —225.8 MPa path. The vertical arrows indicate P< 
and P > (see Fig. [|J . 



from a thermodynamical point of view, the subsequent 
dynamics is completely different in the two cases: af- 
ter P < the system expands, while after P > it contracts. 
From a landscape point of view, P < and P > differ both 
in the depth and the shape of the sampled basin. The 
picture that emerges is that, at P < , the system populates 
local minima which are systematically characterized by 
energy higher and basins of attraction steeper than the 
ones explored at P>. 

It is particularly important to note that, after about 
20 — 30 ps, corresponding to the time requested to bring 
T and P to equilibrium, the Pis value stabilizes, while 
both depth and shape of the explored basins contin- 
ues to change with time. This suggests that the vi- 
brational component to the pressure (defined as P v u, = 
P — Pis 0, 0|) is rather insensitive to the basin prop- 
erties, in agreement with similar findings from equilib- 
rium studies ^]|. The fact that Pjg is constant offers 
an unique possibility to estimate how and if the aging 
dynamics proceeds via a path which is usually explored 
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in equilibrium, by comparing properties in equilibrium 
along a constant Pis path and the aging dynamics at 
the same constant Pjs value. Indeed, in equilibrium, all 
landscape properties (i.e., eis, Pis, S) are function only 
of T and V. By eliminating T in favor of Pis, &is can 
be expressed as a function of V and Pis- Therefore, 
for a fixed value of Pis, &is becomes, in equilibrium, an 
unique function of V. This equilibrium relation can be 
compared with the same relation along the aging path 
to estimate which part of the aging dynamics follows the 
equilibrium properties and which part is instead different 
from it. A similar comparison can be performed between 
basin shape S and V at constant Pis- 

This comparison is shown in Fig. 0] for the two cases 
where a monotonic U-relaxation is observed (large Ti 
values, T t = 220 K and T x = 280 K), and for the case 
Ti = 1501^ where a maximum in the time evolution of V 
is clearly detected. The important information standing 
out from the comparison is that, while equilibration with 
large T\ values follows equilibrium paths, in the other 
case only the final part of the dynamics follows a path 
which goes via a sequence of states which are explored 
in equilibrium. This confirms that only one (P>) of the 
two states with the same T, P, and V can be related, for 
example via a fictive temperature 0] , to an equilibrium 
state, while the other one (P<) has a structure which is 
never explored in equilibrium. Therefore the two states, 
although having the same T, P, and V, are quantita- 
tively different. In other words, while at P> the system 
samples landscape properties which are sampled by the 
liquid at equilibrium at a higher temperature (which can 
be used as additional parameter to quantify the glass 
properties), at P < the system explores a region of the 
landscape which is never explored in equilibrium, since 



the relation between the landscape properties character- 
izing P < are never encountered in equilibrium. Under 
these conditions, it is not possible to exactly relate the 
glass structure to an equilibrium structure and hence de- 
fine a fictive temperature for the system. 

The present numerical study shows that, only when 
the change of external parameters is small, or when the 
system is close to equilibrium, the evolution of the equili- 
brating system proceeds along a sequence of states which 
are explored in equilibrium. Under these circumstances, 
the location of the aging system can be traced back to 
an equivalent equilibrium state, and a fictive tempera- 
ture can be defined. In this approximation, a thermody- 
namic description of the aging system based on one ad- 
ditional parameter can be provided. When the external 
perturbation is significant, like in hyper- quenching ex- 
periments ll'.j . then the aging dynamics propagates the 
system along a path which is never explored in equilib- 
rium. In this case it becomes impossible to associate 
the aging system to a corresponding liquid configura- 
tion. It is a challenge for future studies to find out if a 
thermodynamic description can be recovered decompos- 
ing the aging system in a collection of sub-states, each 
of them associated to a different fictive temperature — a 
picture somehow encoded in the phenomenological ap- 
proaches of Tool and co-workers [2£| and Kovacs and co- 
workers |2l| — or if the glass, produced under extreme 
perturbations, freezes in some highly stressed configura- 
tion which can never be associated to a liquid state. 
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